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Abstract 

The differential decay spectrum dV/dEedq'^ for the semileptonic decay of an unpolar- 
ized hadron containing a b-quark, and the decay spectrum dV / dEedq^dcos 9 for polarized 
Ab decay are computed to second order in the l/m^ expansion. Most of the I /ml correc- 
tions have a simple physical interpretation, which is discussed in detail. The implications 
of the results for the determination of V^b are discussed. The decay spectra for semilep- 
tonic decay of hadrons containing a c-quark are also given. There is a subtlety in the use 
of the equations of motion at order 1/m^ in the heavy quark expansion which is explained. 
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1. Introduction 



The inclusive lepton spectrum from semileptonic B decays has undergone intensive 
experimental and theoretical study. An understanding of it provides information on the 
weak mixing angles Vet and V^b- Recently there has been considerable progress in under- 
standing the theory of these decays. Chay, Georgi and Grinstein [Q showed that inclusive 
semileptonic B XeVe decay can be treated in a fashion similar to deep inelastic scatter- 
ing. Using a two step process that consisted of first using an operator product expansion 
and then a transition to the heavy quark effective theory, they showed that dV / dq^dEf, 
[q^ = (pe + Pu)'^) when suitably averaged over is calculable. Their leading order result 
agrees with the free b-quark decay picture. The full power of their method becomes ap- 
parent when corrections to the leading order result are discussed. These corrections are 
of two types, perturbative as{mb) corrections and nonperturbative corrections suppressed 
by powers of rrib- Chay et al. pointed out that there are no nonperturbative corrections 
of order Aqcd/?t^6- In this paper, we will compute the h^QY)/m1 corrections. The order 
as(m{,) corrections have been computed previously [@]@||1][|5|. The results of this paper 
can be combined with refs. — to give the inclusive semileptonic decay including all 
corrections to order Aqq^/tb^ and asirrib)- 

At leading order in ^^(mfe) and Aqcd/"^& the differential B X^eVe semileptonic 
decay rate is 



dT _ \V^t\^Glml 



dq^dEe 87r3 



2 



2E, q 



1 + 



e{2E,mi,-q^). (1.1) 



In the approach of Chay et al. the b-quark mass appearing in eq. ( |1 . Ij ) has a precise 
meaning that is provided by the heavy quark effective theory. This is necessary for the 
statement that there are no Aqcd/ mi, corrections to have content. In the heavy quark 
effective theory the strong interactions of a bottom quark with four velocity v are given 
by the Lagrange density 

C = by{iv D)b^ + .... (1.2) 

In eq. (|TT^) the bottom quark field by satisfies the constraint = by and the ellipsis 
denote terms suppressed by powers of l/m^. The relationship between the b-quark field 
in the effective theory and the "full QCD" b-quark field is 

(1 + 1^) • 

by = ^ ^^ €^"^""•^6+..., (1.3) 



where the eUipsis again denote terms suppressed by powers of l/m^. It is the same mj, 
that appears in eq. ( |1.3| ) that is used in eq. ( [Til) . With the heavy quark effective theory 
given by eq. (|TT^) {i.e. no mass term for by) rub is a physical quantity that can, at least in 
principle, be determined experimentally. For example, the form factors for the exclusive 
decay — > KcCUf, depend on A = Ma^ — rrib and the determination of A from a detailed 
study of this decay together with the measured A^ mass gives the b-quark mass to be used 
in eq. (1.1). 

Bigi, Shifman, Uraltsev and Vainshtein P] have performed an analysis of the K^Qj^/ml 
nonperturbative corrections to the lepton energy spectrum dV/dEe in semileptonic B- 
meson decay. Bigi et al. found that these corrections are determined by the two local 
matrix elements < B{v)\by{iDYhy\B{v) > /2ml < B{v)\gbya^"'Gfj_„by\B{v) > /4ml. 
The later is fixed by the measured value of the B* — B mass difference. 

In this paper we extend the results of Bigi et al. and Chay et a/., and compute the 
Aqqy)/ ml corrections to the fully differential decay distribution dV / dq^dE^, for an unpo- 
larized hadron containing a b-quark. We also consider inclusive polarized semileptonic 
decay for the special case of the A;,. A^'s produced in decays are expected to be 
polarized . The differential decay distribution for a polarized A^ has the form 

A{E,,q^)+B{E,,q^)cose, (1.4) 



dq^dE^dcos 



where 9 is the angle between the electron direction and the A;, spin in the rest frame of the 
A;,. For reasons similar to those given by Chay et al. in the case of S-meson decay, there are 
no Aqcd/w6 nonperturbative corrections to the differential decay distribution in eq. (|1.4|). 
The AqQj-j/m^ corrections to A{Ee, q^) are similar to those computed for dV /dq^dE^ in B- 
meson decay. One simply replaces S-meson matrix elements by Af, matrix elements and sets 
the Ab matrix element of hyga^'^G^jjhy to zero. However, we find that the Aq^j^/m^ cor- 
rections to B{Ee,q^) are not characterized by just < Kh{v , s)\by{iDYhy\Kb{v , s) > /2ml. 
The normalization of B^E^^q^) involves another order AgQ^/m^ correction which arises 
because the heavy quark spin is renormalized at order 1/?ti^. 

We examine the physical interpretation of the Ag^j^/m^ corrections. Most of the 
corrections {i.e. all of those involving the matrix element of by{iD)'^bv and some of those 
involving the matrix element of hyga^^G^j^hy) can be interpreted as arising from the fact 
that in the bound state the b-quark has an effective mass that differs from and an 
effective four velocity that differs from . These differences arise, for example, from the 
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motion of the b-quark in the hadron rest frame. Corrections of this type are similar in 
spirit to those included in models for inclusive semileptonic B-meson decay [|2[ — 0. 

Section 2 contains a discussion of the kinematics relevant for semileptonic Hf, and 
polarized decay. The operator product expansion and the transition to the heavy quark 
effective theory are discussed in Section 3. This section contains a lengthy discussion of 
the computation. Readers not interested in the details are advised to skip this section 
entirely. Section 4 contains a brief discussion on the use of equations of motion in time 
ordered products. It clears up some confusion on this subject that occurred in the previous 
literature. Section 5 gives the differential decay rates for unpolarized iff, semileptonic 
decay and polarized decay. Section 6 is concerned with the physical interpretation of 
the ^qcd/^I corrections. Section 7 gives differential decay rates for unpolarized hadrons 
He containing a c-quark and polarized Ac semileptonic decay. Section 8 discusses the 
prediction for dr{B XueV^) / dq^dE^ near the boundary of the Dalitz plot. This region 
is important for the determination of Vuh- Particular attention is paid to the size of the 
region of E^, that must be averaged over before experimental results can be compared with 
theory. Numerical estimates and plots of the lepton spectrum are given in section 9, and 
concluding remarks are given in Section 10. 



2. Kinematics 

The semileptonic decay of a b-quark is due to the weak hamiltonian density 

Hw = -V,-6 ^ q^^PLh e-t^PLye = -V.h ^ J^, hp.. (2.1) 

where Pl is the left handed projection operator |(1 — 75)- and are the hadronic 
and leptonic currents, respectively. The final quark qj can be either a u- or a c-quark. The 
inclusive differential decay rate for a hadron if^ containing a b-quark to decay semilepton- 
ically, H\y Xu,c ez^e is determined by the hadronic tensor 

Vr;- = {2nf J2 ^PH, -q-px) {H,{v, s)\jf \X) {X\ JJ \H,{v, s)) , (2.2) 

X 

where j = u, c. The spin J hadron state |i?&(f , s)) is normalized to , instead of to the 
usual relativistic normalization of 2Mh,,v^ as this is more convenient for the heavy quark 
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expansion. W^"^ can be expanded in terms of five form factors if one spin-averages over 
the initial state, 

W^"" = -g'^'Wi + v^v''W2 - ie^'^'^^v^qfiWz + q^q^'W^ + {q^v"" + q^'v^) W^, (2.3) 

where v is the velocity of the initial hadron, defined by 

PH, = Mh^. (2.4) 

W\ and W2 have mass dimension —1, Wj, and have mass dimension —2, and VF4 has 
mass dimension —3. (The form factor of ref. [|l| vanishes by time reversal invariance.) 
The form factors are functions of the invariants and g ■ f , and will also depend on the 
initial hadron iff, and the final quark mass m^. The difference between the heavy quark 
mass mfc and the hadron mass Mu^ will be important in our analysis, so we have chosen 
to write the form factors in terms of q rather than the rescaled q = qjmh used in Ref. [jl|.0 
The spin averaged differential semileptonic decay rate is 



dV _ \V,^\' G-p 



2 n2 

[W^q^ + W2 {2E,E, - 1^2) + Wsq^ {E, - E,)] , (2.5) 



dq'^ dEe dE„ 27r3 

where Ee and Ei, are the electron and neutrino energies in the Hi, rest frame, q^ is the 
invariant mass of the lepton pair, and the kinematic variables are to be integrated over the 
region < AE^^Ej^. The terms proportional to q^ or q^ in eq. ( |2.3|) do not contribute to 
the decay rate if one neglects the electron mass. 

The polarized A;, has in addition to the five form- factors in eq. ( |2.3| ), nine spin- 
dependent form factors which are 

W^" = -q-s [-^^^Gi + v^'v''G2 - ie^"'''^v^qpG3 + q^q^G^ 

+ {q^v" + q^v^) Gsj + is^v" + s^v^) Gq + {s^q" + s^q^) Gj (2-6) 

The identity 



^ We will use m to denote quark masses and M to denote hadron masses. 
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has been used to eliminate terms in Wg'^ of the form i{q^e'^'^l^'^VaqpSx — u)) and 

VaQ/sSx — (/U ^ z^)). The form factors 6*4, 6*5, Gj do not contribute to the decay 



rate if the lepton mass is neglected. The differential decay rate is 

dr 



dq'^dEf,dE,ydcos 6 



... + 



\V,,\'G 



47^3 



2n'2 

F n 

— COS X 



2Ef>Ei, 



1 2 



+ Gsq^ [E, - E,]^ {E, + E,- q^ /2E,) (2.7) 



+ Ge - AE,E,) - Gsq^ - Ggq"" {E, - E, + q^/2E,) 



where 9 is the angle between the electron three momentum and the A;, spin vector in the 
Ab's rest frame. The ellipsis in eq. ( p77|) denotes the part independent of cos 6* and is one 
half the expression in eq. ( p.5|) so that integration over cos^ reproduces the unpolarized 
decay rate. 

The form factors in W'^'^ and Wg'^ are given by the discontinuities across a cut of the 
amplitudes T^'^ and Tgf'^, 



{2.i 



It is easy to see that ImT^'^ = —nW^'^ by inserting a complete set of states between 
the currents. T^*^ for the polarized A;, is defined similarly and has nine additional spin- 
dependent form-factors, 



-q ■ s -g^'^Si + v^v''S2 - ie^^'^'^^v^qpSs + q^q'^S^ 
+ {q^v^ -f q^v^") ^5] + [s^V + s^'v^) Se + [s^q" + s^q^) S7 



(2.9) 



The analytic structure of T^'^ (or T^'^) as a function of ■ f for fixed q^ was given in 
ref. in and is show in fig 1. There is a cut from ^/q^ < q ■ v < (^ 2Mh ) i^Ht + 9^ ~ -^j )' 
where Mj is the mass of the lightest hadron containing the final state quark qj. The 
discontinuity across this cut (which will be called the physical cut) gives W'^'^ for 
semileptonic Hb decay. In addition there are cuts along the real axis for Mh,, q ■ v > 
\ {{2MHf, + MjY — (f — -^IrJ corresponding to the physical process eVeHb X (where 
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X contains two b-quarks), and for v ■ q < — a/^ corresponding to the physical process 
e'^UeHb X. There is no cut for —\f(p' < q - v < \/q^^ since {q^Y > Q'^ ^'^y physical 
state, so that, in general, the physical cut is separated from the other two cuts along the 
real axis. (This disagrees with ref. [|1|.) For Mj = and q^ = Mfj^ the two cuts on the 
positive real axis are not separated. The end of the physical cut at v ■ q = Mh^ coincides 
with the beginning of the second cut. Similarly for q^ — the cut on the negative real axis 
ends at the same point the physical cut begins. 

The amplitude T'^'^ can be computed in QCD perturbation theory away from the cuts 
along the real axis. Provided q'^ is not too near zero or M]j^ (for qj = u), the value of the 
amplitude in the physical region can be obtained by performing a contour integral along 
the closed contour C shown in fig. 1 that stays away from the cuts in the complex q ■ v 
plane. 



3. The Operator Product Expansion 

The amplitude T^^'^ can be computed reliably in perturbative QCD in a region that is 
far from the cuts and therefore free of infrared singularities. The time ordered product 

-i j rf^e-^'?-"T(J'^V") (3.1) 

can be computed using an operator product expansion in terms of operators involving 
b-quark fields in the heavy quark effective field theory. The coefficients of the operators 
in the operator product expansion are determined by evaluating the matrix element of 
the time-ordered product between quark and gluon states. Once the operator product 
expansion has been computed, one can compute T^^ by taking the matrix element of the 
operator product expansion between hadron states. We will compute the matrix element 
between unpolarized hadron states to determine the form- factors T^, and we will take the 
matrix element between polarized A;, states to determine the spin-dependent form-factors 
Si. The operator product expansion can be written as an expansion in inverse powers of nib- 
The expansion of T^^ to order 1 / can be written in terms of gauge invariant operators 
of dimension less than or equal to five. The operators will involve the h field, covariant 
derivatives D, and the gluon field strength tensor G^'^ . The coefficients of the operators 
involving h and D are determined by taking quark matrix elements of eq. ( p.l|) , and the 
coefficient of operators involving G^^ are determined by taking gluon matrix elements of 
eq. ( |3.1| ). We only calculate the form-factors Ti_3, S'i_3, 5*6, 5*8 and 5*9 which are the only 
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ones that contribute to semileptonic b-decay when the mass of the lepton in the final state 
is neglected. 

The quark matrix element of eq. between b-quark states with momentum rribV + k 

is 

2—- u Pl {m^i -i + lfe + m,) 7" Pl (3.2) 

[nibV — q + k)^ — rrij + le 

from fig. 2(a), where u is the quark spinor. The crossed diagram of fig. 2(b) has no 
singularities inside the contour C of fig. 1, and does not contribute to semileptonic b- 
decay. The matrix element eq. ( |3.2| ) can be simplified using the identity 

The l/m^ expansion is given by expanding eq. ( |3.2|) in a power series. The momentum q 
can be of order mf,, but k is only of order Aqcd- Thus each factor of k in the expansion of 
eq. ( |3.2| ) corresponds to a l/m^ suppression. The factors of k in the matrix element will 
become factors of iD in the operator product expansion, so each factor of iD corresponds 
to a 1/mf, suppression. 

3.1. The /c° Terms 

The order k^ term in the expansion of eq. (pi. 2]) is 



I {jUhV - q)^ Y + {"f^bV - qY Y - {^bi) - i) g^" 
i\q y. 

-^e^^"^(mb^;-Q)^7/3 ]Plu, 



(3.3) 



where 

Aq = {rribV — qY — + ie. (3.4) 

The matrix elements of the operators b'j'^b and 67^75 b between b-quark states are uj'^ u 
and U'j^'y^u respectively, so the operator product expansion is obtained by replacing u 
and u in eq. ( p.3| ) by the fields b and b respectively, 

^lirribV - qf g"^ + {m^v - qf g^^ - {mbV - qf g^" 



Ao 



(3.5) 

ie^-"^{mhV-q)^j b^xPLb. 



The spin averaged matrix element J2s i^biv, s) \ 67'^75 b \Hi,{v, s)) between hadrons Hh of 
velocity v is zero. The matrix element 

{H,{v,s)\b'yH\Ht{v,s)) = v\ 

to all orders in l/mj,, since b-quark number is an exact symmetry of QCD. (This relation is 
true because we have normalized our hadron states to v^, and defined v in eq. ( |2^ ) to be the 
velocity of the hadron Hb.) Taking the matrix element of the operator product expansion 
between spin-averaged Hb states, and comparing with eq. ( p.8|) gives the contribution to 
T^'^ of terms in the expansion of eq. ( p.2| ), 

Ti'^ = rUb, (3.6) 

^(0) _ _J_ 

' "2Ao- 

The terms in the operator product expansion give b-quark operators that have zero 
derivatives. In principle, these operators could produce contributions to T'^'^ of higher 
order in l/rub because of l/m^ corrections to the matrix element of the operator between 
hadron states. However, the only matrix element we need is the zero momentum transfer 
matrix element of a conserved current, which has no 1/mf, corrections. 

The matrix element of the operator 67^756 between polarized A;, states does not 
vanish. At leading order (in l/nib) it is equal to the spin vector s^. However, the axial 



current is a generator of heavy quark symmetry [|T0[ that is broken by 1/ rub terms in the 
effective Lagrangian. Thus the matrix element of the axial current between polarized A;, 
states is not equal to the spin four vector. The first correction to the matrix element is of 
second order in symmetry breaking [|11|, so we write 

{Ab{v,s)\b'j^'j5b\Ab{v,s)) = (l + eb) u{v, s)'j^'j5u{v, s), 

(3.7) 

= il + eb)s\ 

where eb is a parameter which is of order Aq^j-j/m^, and is defined by eq. (|3.7| ). Substituting 
this into eq. (|3.5|) gives the spin-dependent form-factors from the k^ terms 

(3.8) 

o(0) _ '^(, ^(0) _ l.-i , N 1 

si''> = ) = 0. 
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3.2. The Order k} Terms: Spin Averaged Case 
The linear terms in k in eq. (|3.2|), 



2k ■ {rribV - q) _ , f i i a\ (3 9) 
—2 u<[mhV — q) '^ + [mbV — q) 'J — [nibf — q) g ^ 

Aq I 



produce operators in the operator product expansion with one derivative, of the form 
h"^^ iD'^ b and 67^75 iD'^ b. The matrix elements of these operators need to be computed 
to first order in l/rub, since they contribute to terms that are already suppressed by 
l/rub. Unlike the current operator b'~^'^b in the k^ terms, these operators will have l/rub 
corrections to their matrix elements, and so can contribute terms to T^'^ at order l/m^ 
and l/m^. The matrix element of 67^75 iD'^ b vanishes between spin averaged Hb states. 
Its contribution to the polarized Af, decay amplitude is discussed in the next subsection. 
The matrix element of b'^'^iD'^ b can be computed in a l/rub expansion using the heavy 
quark effective theory. 

The b-quark is represented by the velocity dependent b-quark field by in the heavy 
quark effective theory. Since the terms linear in k are already of order l/m^, we only need 
the relation between b{x) and 6„(x) to first order in l/rub 



b{x) = e~ 



■imf,v-x 



iV) 
1+ ^ 



by{x). (3.10) 



2mb 

The QCD lagrangian for the b-quark in the heavy quark effective theory is 



C = byivDby + by by - Zb by b^ + o(^]. (3. 11) 



where Zb is a renormalization factor, with Zb{fi = nib) = 1- The operator by [iD)^ by/2nib 
is not renormalized because of reparameterization invariance |T^ . 



Eq. (|3.10| ) gives the expansion of the operator in the effective theory 



67^ iD-'b = by 7^ iD-'by + by ^^^7^ K + by 'J^ iD'' ^ by, 

2nib 2nib 

= v% iD^by + —by iD^^ iD^^ by —by g G'^^ o by, 

nib 
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(3.12) 



where we have used the relation by ■y'^^by = v^by by (which foUows from the constraint 
= by) and the commutator [D^,D'^] = igG^'^ . In eq. (3.12) the brackets around 

indices denote that they are symmetrized. 

The spin averaged matrix element of the first term of eq. (|3.12| ) must have the form 



(2J + 



f , s) \ by iD'^ by \Hi,{v, s)) = Av'^ , 



(3.13) 



where A is a constant to be determined. Contracting both sides of eq. ( p.l3| ) with f ^ gives 



A = {Hb{v,s)\byiD ■ vby \Hi,{v,s)) 



(3.14) 



The coefficient A is zero at lowest order in 1/mb, since {v ■ D) by = is the lowest order 
equation of motion in the heavy quark effective theory. Thus the order k terms make no 
contribution to T^^'^ at order l/mf,. The order terms also did not contribute to T^^'^ at 
order 1/mb, so there is no 1/mb correction to T^^'^ . This reproduces the result of ref. 
The parameter A is non-zero at first order in l/rub, 



A = -{Hb{v,s) 



bi, ^ — ^— by — Zb by — — br. 



2mb 



Anib 



\Hb{v,s)) 



(3.15) 



using the equations of motion from the Lagrangian eq. ( 3.11|) to order l/rrib- This l/nib 
value of A contributes to T^*^ at order 1 / . 

It is useful to define the dimensionless parameters 



Eb = -(Hb(v,s) 



Oy — — Oy — Zjb Oy ; 0, 



2ml 



Ami 



\Hb{v,s)) 



Kb = -{Hb{v,s)\by 



2ml 



by \Hb{v, s)) 



(3.16) 



Gb = Zb {Hb{v, s)\by ^^"^Z — by \Hb{v, s)) 

Ami 



to characterize the 1/m^ corrections, with Eb = Kb + Gb- Gb, Kb and Eb can be thought 
of as the average value of the spin-energy, the kinetic energy, and the total energy of the 
b-quark in the hadron Hb, in units of mb- AH three parameters are expected to be order 
AqQ^/m^. The operators in eq. ( |3.16| ) are renormalization point independent, since the fi 
dependence of Zb cancels the n dependence of by gG^fBcr"'^ by. The anomalous dimensions 
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of the operators only affect the relations between the parameters for different heavy quarks. 
For example, the parameters for b and c quarks are related by 

mlKc = mlKb, 

(3.17) 

mlGc/Zc = mlGb/Zb. 

Since Zcijn^ = 1 and Z\y{mb) = 1, the ratio Z^jZ^ is given by the scaling of 6„ gGapcr"'^ by 
between the scales rric and m^, Z^jZ^ = \as{jnh)las{vn^^^'^'^ [0. The gluon operator G^" 



in the operator product expansion occurs through the equations of motion, and through 
the commutator ^ = igG^'^ . We will use to parameterize the matrix elements of 
the gluon operators obtained using the equations of motion, and Gb to parameterize the 
matrix elements of the gluon operators obtained from [D^,D'^]. This distinction will be 
useful in sec. 6. With this convention, 

A = rubEb, (3.18) 
and the matrix element of the first term of eq. ( pj.l3D can now be written as 



-^^^ J2 ^) I ^)) = ^bEbv\\ (3.19) 



The matrix element of the second term in eq. (|3.12|) must have the form 

j^^^^ J2 ^) I '^^^ '^^^ l^''^^' ^)) = (^1^^^" + B2v\^) . (3.20) 

This term has two covariant derivatives, so we only need its matrix element to lowest order 
in l/rub- The lowest order equation of motion [v ■ D)by = implies that Bi + B2 = 0. 
Taking the trace and comparing with eq. (|3.16|) gives 



s 



The operator by gGaT(^°'^ by in eq. ( p.l2| ) is renormalized at a scale n = nib, since 
that is the scale at which the operator product expansion has been performed. We can 
therefore multiply the operator by the renormalization factor Zb, since Zb{mb) = 1. This 
makes the operator renormalization group invariant, and includes the QCD scaling of the 
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operator between m^, and fi. The matrix element of the third term in eq. ( |3.12|) must have 
the form 



2J + 



^ J] {Htiv, s) I K gG'^'^aJ 6„ \H,iv, s)) = {C,g^^ + C2v\'') . (3.22) 



Contracting both sides with v'^ , and using by a'^'^Vr 6„ = gives Ci + C2 — 0. The trace 
gives 



2J + 



1 



(3.23) 



on comparing with eq. ( |3.16| ). 

Substituting eq. ( |3.19| ), ( |3.21|) and ( |3.23| ) into eq. (|3.12|) , and substituting the result 
into eq. (|3.y|) gives 



1 {mb — q - v)' 



2Ao 



IS? 
^0 



+ ^iKb + Gb) -^ + - 



2 (q-v)^ 



^0 



Ao A? 



^0 



+ ^ (i^, + Gb) 



1 2mbq ■ V 



(3.24) 



(1) 



-TUbEb 



{nib - q-v) 



^0 



2mb 



{Kb + Gk 



nib — q ■ V 

A2 
^0 



3.3. The Order Terms: Polarized Ab Case 

The spin-dependent form-factors arising from the order k terms in the operator prod- 
uct expansion of eq. ( |3.9|) are obtained by taking the matrix element of the operator 
67^^75 zD"^ 6 between polarized A;, states. This operator did not contribute to the spin- 
averaged matrix element between unpolarized Hb states discussed in the previous section. 
The method used to evaluate the matrix element is similar to that used for the operator 
b'j^ iD'^ b in the previous subsection. The operator can be written in terms of the field by 
of the effective theory, 

67^75 iD-" b = by -f^-f5 by + -^by ^ iD^ by + -^by ^f^^f^ by. (3.25) 

znib Irrib 

The matrix element of the first term of eq. ( |3.25| ) between polarized A;, states has to have 
the form 

{Kb{v, s) I by 7^75 iD-" by \Kb{v, s)) = A 117^75 u v\ (3.26) 
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by heavy quark spin-symmetry. Contracting both sides with v'^ and using the equations of 
motion determines A = uibEb, where Ef, is defined in eq. The matrix element of the 

second and third terms in eq. ( p.25|) can be simphfied by neglecting terms proportional to 
QaP_ fpj^g operator G"^ vanishes in any matrix element between A;, states at zero recoil, 
since the light degrees of freedom in the Af, have spin-zero. (That is why spin symmetry can 
be used in eq. ( |3.26| ) even though we are including effects of order l/mt.) The only vector 
that can be constructed using the light degrees of freedom is f^, and it is not possible to 
construct a tensor that is antisymmetric in two indices from a single vector. 

The 1 / nib terms in eq. (|3.25[) can be simplified using 7- matrix algebra and neglecting 
G"^ to give 

'I 

which is equal to 



-Aa)75&., (3.27) 
nib 



-J_b^iD("iD^K^-f^-f5K, (3.28) 
nib 

using the equations of motion. The matrix element of eq. ( |3.28|) between polarized A5 
states is 

-— {Ab{v,s)\KiD^''iD^^v^-f^-f5bv\Mv.s)) = \mbKbV^uYl^u, (3.29) 
nib 3 

using heavy quark spin-symmetry, and the matrix element eq. ( p.21| ). Eqs. ( |3.26| ), ( |3.27| ), 

and ( |3.29| ) imply that one can make the substitution 

- 2 

67^75 iD'' h nibEb s^v^ + -nibKb v^s^, (3.30) 
in eq. ( p.9|) to obtain the spin-dependent form factors 

S[^^ = ^{nib-q-v) (^nibEb + '^nibKb^ , 
= i^rnlKb, 



3Ao 



S, 



(1) 




6 ~ 


2Ao ( 


(1) 

8 ~ 


2Ao (" 


r(l) 

9 - 


1 . 
^(mt 

^0 



(3.31) 



3 ' V 
2 „ \ 1 
^0 



^s' = 777— ( '^bEb - ^nibKbj - -^{nib-q-v) mjEb, 
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3.4- The Order /c^ Terms 



+ 



The order k"^ terms in eq. ( p.2|) are 
2fc-(m^z;-g) _|^^^, ^ ^^^^ _ ^^^^ _ 

H I (m{,t^ — q)'^ Y + (mft-y — qY Y — {trtbi) ~ i) 9^"" (3-32) 



{k ■ {nibV — q)) k'^ 

i-in i-in 



The matrix element of eq. ( |3.32| ) between unpolarized hadrons Hb can be written in terms 
of the operator b'j'^ iD^^iD^^^ h. This matrix element is needed to lowest order in l/m^, so 
h can be replaced by the heavy quark field 6^, and 7^ replaced by f ^. The matrix element 
needed is 



2jTT ^ iD^-'^zD^^ K \Hb{v, s)) = (g-^ - v<^vP) , (3.33) 



s 



using eq. (3.21). Substituting eq. ( p.33|) into eq. ( p.32|) gives the contribution of the k^ 
terms to T'^^ , 

Ti^ = -KnlKb{rnb- q-v)\^-^^ [q^ -{q-vf] -^}' 

The spin-dependent contribution to T^^ for polarized A5 states is given in terms of 
the the matrix element 

(A,(^, s)\K 7^5 iD^'^iD^K \^h{v, s)) = -^mlKb (^7°^ - v'^v") s\ (3.35) 

using eq. (|3.21| ) and heavy quark spin-symmetry. This gives the order /c^ contribution to 
the spin-dependent form-factors 
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^0 



= = 0. 

3.5. The one gluon matrix element 

Finally, one needs to compute the one gluon matrix element of eq. given in fig. 3 
to determine the coefficient of the G"^ operators. The one gluon matrix element can be 
expanded in a power series in the momentum p of the external gluon. The terms of order 
are identical to the one gluon matrix element of the operators we have already found, 
with the gluon field from the covariant derivative D = d -\- igA. The terms linear in p are 
the matrix element of the operator 



9 T 



2AI 



bG'^^e^pxAmbV - q)^ [g^'^Y + ^7""7^ - ^/''"t" + ^e'^"""7r] Pl K (3.37) 



This is a dimension five operator, so we need its matrix element to lowest order in l/m^. 
The fields h can be replaced by the heavy quark fields 6„. The matrix elements needed to 
evaluate the one-gluon contribution to T^'^ for a spin-averaged hadron Hh are 

^ {Hb{v, s) I K gG'^P-i^ K \Ht{v, s)) , (3.38) 



^ J2 {Hb{v, s)\K ^/G-^T^s K \H,{v, s)) . (3.39) 



2J + 

and 

2J + 

The matrix element in eq. (|3.38| ) can be simplified by replacing the 7 matrix between by 
fields by v. The resultant matrix element must vanish, because there is no antisymmetric 
tensor in the indices a and (3 that can be constructed out of the single four vector v. The 
matrix element in eq. (^.39|) must have the form 



^ J2 (Hbiv, s) I by ^7G°^7^75 by \H,{v, s)) = Ne^^^^v-,. (3.40) 



2J + 

s 
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The constant can be evaluated by contracting both sides of eq. ( |3.40|) with ea/sxpV^ to 
give 



2J + 



^ J2 ^"/^V^"" {Hb{v, s) I 6„ A^75 K \Ht{v, s)) = -6N. (3.41) 



s 



Using the spinor identity 



(^apXpV^ by 7^75 by = -by ac,f3 by, 



and eq. ( |3.16| ) gives 



^^(if,(^,s)|6.5GA'75&.|i?;,(^^,s)) = ^m^G.e-^^X. (3.42) 



2J + 

s 



Substituting eq. ( |3.42D into ( |3.37| ) gives the contribution of the one gluon operator to T^'^ , 

4'^ = l^nlG,^^, (3.43) 
7^3 --g^'^^^Af 



The gluon terms do not contribute to polarized decay since G^^ has vanishing matrix 
element at zero recoil between Af, states. 



3.6. Summary 

The final expressions for T'^'^ and T^^ to order l/mf are obtained by combining 
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eqs. (H), O, (EH), (pTI), (PI), (PI), and (P|) 



To 



5i 

^3 

^8 
Sq 

where 



1 2 

—-{mb~q-v){l+ Xb) + -rub {Kb + Gb) 



2A 



1 ^f-{Q- v) 



^0 



^(l + X,) + |m,(K, + G',) 



1 2m;,q' ■ V 



+ 



TUbEb 



+ o^bKb^ + i^rn^Gb^ 



4 

-■; 
3 







^ (1 + Xb) - \mb {Kb + Gb) ^^^^^ 



1 1 



'^iGb-To, 



1 

A3 

^0 



3A2 



3A2 



A3 

^0 



A3 

^0 



-mlKb[q^-{q-v)^] 



hmbKbq-v 1 1 
2(^ + ^^^A^' 



ml Kb, 



rubKb 



-m 



\KbW-{q-v)'] 



^m%Kbq-v brUbKb 1 .rUb 



3Ag 



-m 



]Kb[q'-{q-v)'] 



+ 



bmlKbq ■ v rubKb 1 



rub 



3Ag 



A3 

^0 



-mlKb [q^ - {q ■ v)^] 



Kb 
A2 

^0 



"mbq ■ V + -m^ 



Xb 



mbEb 



8 oK 



{rub-q-v) -^ml—^ q^ - {q ■ v 



3 "Al 



1. ^ 1 

2^^ + ^^)a^' 



(3.44) 



(3.45) 



The expressions in eq. ( p.44|) can be simphfied using the identity Eb = Kb + Gb- There 
is an additional simphfication in A;, decay, where Gb = 0. 



4. Time Ordered Products and the Equations of Motion 

In the above computation, we used the equation of motion to order 1/mb. An- 
other method commonly used in the literature is to treat the l/rnt terms in the La- 
grangian eq. (|3.11| ) as a perturbation, so that the equations of motion is {v ■ D)by = 0. 
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The 1/mi, terms in the Lagrangian then give terms that are time-ordered products 
of operators with 1/mi, terms in the Lagrangian. For example the matrix element 
{Hh{y,s)\byiD'^hy\Hh{y,s)) (that was needed in sec. 3, see eq. (|3.19|) ) is zero using 
{v ■ D)by = 0. However, one now gets an additional contribution 

{Hb{v,s)\TS^KiD^b,ix) i J A^^l \Hb{v,s)), (4.1) 

where the l/m^ terms in the Lagrangian are denoted by Ci/2mi,, using the notation 
of ref. The matrix element eq. (|4.1| ) is equal to Av'^ , where A is the constant of 

proportionality. Contracting both sides with v'^ gives 

A = {H,iv, s)\T I^Kiiv ■ D)K{x) I j A^^} \Hb{v, s)) . (4.2) 

C-i{y) is a sum of operators of the form bv{y)Xby{y), so that the time ordered product in 
eq. (U) is 



J dy T {by{iv D)by{x) byXby{y)] = j dy T {by{iv D)iS{x,y)Xby{y)} , (4.3) 

where S{x,y) is the b-quark propagator, which satisfies the Green function equation {iv ■ 
D)S{x, y) = d{x, y). This converts the T-product into the local operator byXby, so that A 
is given by 

A = -{H,{v,s)\^\H,{v,s)), (4.4) 
the same answer as that obtained by using the equations of motion to order 1 / rrib . 



5. Decay Distributions for Hadrons Containing a b-Quark 

The amplitude W^'^ can be determined by performing a contour integral of T^^ from 
eq. ( |3.44|) around the contour C of fig. 1. This contour integral is trivial to do, and is 
equivalent to taking the imaginary part in T^'^ directly by making the replacements 



1 

1 

^0 



A3 2 
^0 



6 [{mbV — — m^) 

-5' {{rribV - qf - m^) (5.1) 

^5" {{rrn,v - qf - m") 
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in eq. (|3.44 ). This is very different from the analogous calculation in deep inelastic scatter- 
ing. In our problem, there are non-trivial function of q and f , multiplying a few operators, 
so that the amplitude has the form fi{q^ , q-v) {O)-. The QCD corrections multiply each op- 
erator by an anomalous dimension, fi{q'^, q ■ v) (O)^ \ifi{q^ , q ■ v) {O)-, so the imaginary 
part of T^'^ is related directly to the imaginary part of fi. In deep inelastic scattering, the 
different powers of g ■ f in the expansion of fi would each be multiplied by the anomalous 
dimension of a different twist two operator, so that the imaginary part of T'^'^ is related 
to the imaginary part of fi by a non-trivial convolution. 

One can now compute the inclusive lepton spectrum by substituting eq. ( |3.44| ) and 
eq. ( |5.1| ) into eq. ( |2.5| ), and integrating over Ei,. The kinematic region is determined by 
the mass of the decaying hadron Mh^- The 5 functions in Wi restrict the integration to 
that given by the parton model kinematics determined by the quark mass m^. In QCD, 
one can prove that M//,, — 'rrib = A > 0, so that the parton model kinematic region 
is contained within the hadron kinematic region. Thus there is no dependence on the 
hadron mass Mh^, through the limits on the region of integration, and the lepton spectrum 
is determined only in terms of the quark mass m^. This would not he the case if A. were 
negative. Evaluating the E^, integral gives the decay distribution for Hb Xez^e, 

^ '^^ Q{z)[\2{y~e){^ + e-p-V) 



(5.2) 



where 



Fb dy dq^ 

+ l2Eb{2q' - 2fp + y- 2fy + py) + SKb{2q^ - (f + fp - 3y) 
+ 8Gt{-f + 2ct - 2fp -2y- 2fy + py) 

+ 8{z)Wl2Ebq\y - + 2y- y^) 

y I 

+ AKb{-cf + 9qS - Qfy^ - 2qV " fv^ + v') 

+ 8fGb{y - + 2/ + + Kb5\z)'^{y^ - ffiy - f) 

J 



1 + f-p-^-y, (5.3) 



^2 ™j 2£;e 

q = p = ^, y = — , 5.4 

ruh 



and 
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+ EAA{?,-y)y' + 



2 , 12yV 4(6 - 4y + 2/2)^2^3 



The curve 2; = is an edge of the the Dahtz plot in the (f' — y plane. Integrating eq. ( 5.2 
with respect to gives the inclusive lepton spectrum, 



(5.6) 



3 (l-y)2 {1-yY 3(1 -y^ 

This agrees with the result of Bigi et al. Integrating eq. ( ^T6|) with respect to y gives the 
total decay rate 

^r= |l-8p + 8p3-p4_i2p2iog^ 



4y^9 + 2y) , 4y2(2|/2 - 2y - 3)^^ , 4y^ {18 - lOy + 5y^ - y^) 

^""^^ ^ — ^ 0^^? ^ 

4y2(15 + 2y) , 8yH3-2y)p , 12yV , 8?/2(-6 + 4y - y^)^^ 
+ (jrbS 1 7z r?; r — + 



+ Eb^5-24p + 24p2 _ 8p3 + 3p4 _ ^2^2 log p 
+ Kb I -6 + 32p - 24p2 _ + 24p2 log p 
+ G J -2 + 16p - 16p2 + 2p4 + 24p2 log p 



(5.7) 



which also agrees with Bigi et al. 

The decay rate for 6 ^ -u is given by eqs. ( |5.2| )-( |577|) in the limit that p ^ 0. This limit 
must be taken carefully because of the 1/(1 — y) singularities in eq. ([5.6|). The resulting 
expressions are 

^ 9{z)ll2{y - g2)(i + f _ ^) + 12£;,(2g4 + y - 2fy) 



Tb dy dq^ 

+ 8^5(2^"' - t - 32/) + 8Gb{-f + 2q^ -2y- 2fy) 

+ 6{z)\\l2E,f{y - + 2y- y^) (^.8) 

+ 4X,(-Q^6 + ^cfy - Qfy^ - 2qS^ - fy^ + y'>) 
+ SfCiy - + y + 2/')} + d\z)K,^{y^ - ffiy - g^^), 

where now ^ 

^ = 1 + y, (5.9) 

y 
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and 



1 dr 

To dy 



2(3 - 2y)y^ + 4(3 - y)y^ + 25(1 - y) 



^Kb^^'-y^-lsil-y)^ls'il-y) 



(5.10) 



+ Gt 



4|/2(i5 + 2y) 16 



+ -S{l-y) 



r = r, (1 + 5Eb - 6Kb - 2Gb) . 



(5.11) 



The results eqs. ( p.lO|) and ( |5.11|) agree with Bigi et al. The (^'-function in eq. (|5.10|) is 
present because the parton model decay distribution oc 2(3 — 2y)y^ does not vanish at the 
end point. This will be explained in more detail in the next section. 

The decay distributions for a polarized A;, have the form A + B cos 9. The coefficient 
A is half the value of the corresponding decay distribution for an unpolarized given in 
eqs. ( |5.2D -( |5TTT|) . The results can be simplified for A^ by setting Eb = Kb and Gb = 0. 
The coefficients of the cos 9 terms are 



1 



dV 



Tb dy d(p dcos 9 



. . . + 



9{z) 



6 (1 + eb) {y - f){-2q' + y + fy - py - y'') 



+ 2Kb{-Qct + I2^'y + ^''y - 4q''py - 3y' - 6fy^ + 3py 



+ 5{z)^{Kb{-Af + 2qS + 3qV - qV - + r 



Kb5'{z)%iy'-qY{y-f) 



COS 9, 



(5.12) 



1 



dV 



Tb dy dcos 9 



+ 



(1 + e,) (l-2y)y2-3yV + 



3y'p' {l + y)y^p^ 



+ Kb 



10y3 ^2y3(5-2y)p2 4y^{5 + 2y - y^)p^ 



3(l-2/)5 



cos 9, 



(5.13) 



1 dV 

Th dcos 9 



. + (1 + ef, - Kb) 



i+2p + 6p^-fp3_1^4 



+ 6p log p + 4p log p 



cos 9. 



(5.14) 
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The p — s> limits of the polarized Af, distributions are 



1 



dV 



T(, dy d(p dcos 9 



0{z) 



6 (1 + £5) {y - f){-2f + y + fy - y^) 



+ 2irfe(-6g* + I2qi'y + Aq'y - 3y' - 6fy^) 



+ 5{z) - { K,{-4e + 2fy + s^y' - q'y' - q'y^ + yl 

COS 9, 



y 

K,S'{z)%{y'-my~Q') 



1 dV 
Tb dy dcos 9 



1 + e,) (l-22;)r 



1 dr 



1 



cos 6*, 



{1 + eb- Kb)cos9. 



(5.15) 



(5.16) 



Tb dcos 9 6 

The formulae obtained here will be discussed in detail in the following sections. 



(5.17) 



6. Physical Interpretation of the 1/m^ Corrections 

There is a way to obtain most of the 1 / corrections which also provides a physical 
picture of these corrections. One can obtain a class of 1/m^ corrections by taking the 
lowest order expression eq. ( ^.6) ) for T^'^ (or W'^'^) and smearing it over a distribution of 
b-quark momenta in the Hb hadron. This will give the terms proportional to Eb and Kb 
in sees. 4-5, but not those proportional to Gb- This is why we did not simplify the results 
using Eb = Gb + Kb. 

We begin with the spin independent terms. The b-quark momentum in \Hb{v, s)) can 
be written as p = rubV + k. The lowest order (parton model) expressions Tq'^ {q, v, nib, rrij) 
were obtained by considering the decay of a quark of mass rub and velocity v in the rest 
frame v = (1, 0, 0, 0). A quark with momentum nibV + k can be considered to be an on-shell 
quark with mass and velocity v' , where f'^ = 1, and m^f' = rubV + k. The decay rate 
of such a quark can be obtained by using Tq'^ [q, v' , m^, mj)/v'^, where v'^ is the Lorentz 
time dilation factor for a moving particle. The value of T^^'^ is obtained by averaging 
TQ^{q, f m^, mj)/v'^ over a distribution of b-quark momenta. This average is most easily 
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done by writing f'^ = v ■ m'^v' /m'^. m'^TQ'^ {q,v' ^m'^^nij) /v ■ m'^^v' can be written as a 
function of q, rrij and the product m'l^v' . This makes the averaging simple, because one 
can use the substitution m'fjV' = mi,v + A; to rewrite the expression in terms of rribV and /c, 
and then averaging over k. Since q is unaffected by the averaging, terms proportional to 
qt^ or q"^ will remain proportional to q'^ and q'^. Thus averaging T4 or T5 will not produce 
Ti_3. Similarly, averaging Ti or T3 will produces terms only proportional to Ti or T3, but 
averaging T2 produces terms proportional to both T2 and Ti. 

As a simple example, we will consider the average of the piece of T^*^ containing T3 
explicitly. The average we need is 



<_^.^,f^'^^P^'^q^Ts{q,q-v',mi)) 



. I / (6.1) 




V ■ {rribV + k) \ {rribV + k — q)"^ — m 

The averages over k can be written in terms of 

{k^kf) = \ {k") ((7"^ - v'^vf') = -iK^ml [g^P - v'^v^) 




(6.2) 



where Ei, and K}, are the mean total energy and mean kinetic energy in units of m^. This 
gives the terms in eq. (|3.44|) with the exception of the Gh term. A similar computation 



also reproduces Ti and T2 in eq. ( p.44|) . 



The above averaging procedure computes the 1 / corrections using a distribution of 
quark momenta /c in a hadron Hi. Here k is considered to be a number, not an operator. 
The operator product expansion gives a similar result, with k replaced by the covariant 
derivative iD. The commutator of two covariant derivatives is proportional to the gluon 
field-strength, [D",!)^] = igG""^. The covariant derivatives commute if we neglect all 
terms involving G"'^, i.e. all terms involving Of,. Thus the 1/mf, terms obtained by the 
averaging method are identical to those obtained using the operator product expansion 
neglecting G"^. This result holds to all orders in the l/rrib expansion due to reparameter- 



ization invariance 13 . 
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Our averaging procedure provides a useful check on the total decay rate. The total 
decay rate (neglecting 1/mf, corrections) can be written as Fq = m^/(m^/m^). The total 
decay rate for a quark distribution is given by averaging 

rriij j(m jm^ )/v = j{m jm^ )/v ■ rriijV . 

Now m[,2 = m'f^'^v'^ = {rubV + k)^ = m^(l + 2Eb - 2Kb), and v ■ m[v' = mb{l + Eb). This 
gives 

r = ml{l + + 2Eb - 2Kbff{p{l + 2Eb - 2Kb)-'), (6.3) 

where p = m? / . Expanding this and retaining corrections to order 1 / gives 

r = To + Eb{hTo - 2p^) + Kb{-%T^ + 2p^). (6.4) 

dp dp 

The result eq. ( |5.7| ) agrees with this check. 

One can also use the averaging procedure to determine the leptonic spectrum. Recall 
that we needed to compute the average of tti'i^Tq'^ /v -m'^v' , where m'^TQ^ is a function only 
of m^f ' and q, and has mass dimension zero. This can be written as 

V ■ m',v' / ~ mb{l + Eb)Y'^ ^ dmbV^^'^ 
1 

+ fc" hf rribT'"' 

2 dmbV^drubvl^ ^ 

Tr-EbTr + EbV-^rubTr 



(6.5) 



1 

- -KbTHb (g""^ - v'^v^) - mbT^" 

3 drribV^'dmbvl^ ° 

using eq. ( |6.2|) . The decay rate depends on L^i,T^'^ = F which has mass dimension two, 
and is a function only of kf,, k^,, mv, rrij and q. One can rewrite F in terms of the 
dimensionless variables y, q p, and an overall factor of m^. The variables y, q and p were 
defined in eq. (|5.4|), and x is defined by 

2E„ 2mbV ■ ky 
X = = — —. (6.6) 

rrib ymbv)'^ 

The averaging formula eq. (|6.5| ) for T^'^ implies that F can be written in terms of the lowest 
order expression Fq = L^j^Tq"^, by differentiating with respect to rubV. In performing the 
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differentiation, it is important to remember that Tq^ and Fq are to be considered as 
functions only of the product m^f , not of mj, and v separately. Thus nib is an implicit 
function of f , with 



d 



d 



Using the partial derivatives 
dx 



-nib 



dnibV^ 



[{mbV) ■ {nibV)]''^'^ = Vol- 



(6.7) 



d{mbv'^ 

dy 

d{mbV'^ 



ml 



nib 



mf 



rrib 



d{nibv'^ 

dp 
d{nibV'^ 



2^ 

Vc 

rrib 
nib 



(6.8) 



and eq. ( |0| ) our averaging procedure implies that 



l + Eb 
+ Kb 



l-2p—-2f — 



d 

- — X- 



d 



dp d(p dy dx- 



2 d 

2 + 2p— +2f— + 2y— + 2x 
dq^ 



d_ 
dp 



d_ 

dy 



d_ 

dx 

C>2 



(6.9) 



^ \y2^ ^ ix^-^ + -(xy - 2q^ 

3 dy'^ 3 dx'^ 3 dxdy 



Fn. 



The differential decay rate is given by integrating F over the phase space, which is pro- 
portional to dydxdq^. Integrating eq. (|6.9| ) with respect to x gives the formula 



dV 



dy d(p- 



l + Eb 
+ Kb 



2 2^ 2^2 ^ ^1 
dp dcp dy 

— + 2 — + 2q^— + -y— + -y'^ — 
3 dp dq^ dy 3 dy'^ 



dVr 



(6.10) 



dy d(p 



where dTo/dydq'^ is the parton model decay rate obtained by setting Eb^Kb^Gb — > in 
eq. ( p.2| ). Integrating eq. ( |6.10| ) with respect to q^ gives 



dr 

dy 



l + Eb 
+ Kb 



4 _ 2p— - y — 
dp dy 
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d Ad I ^ d"' 



3 ~^ '^^ dp ^ 3^ dy ^ 3^ dy'^ 



dTo 
dy 



(6.11) 



Integrating eq. (|6.11|) with respect to y gives 

d 



l + Eb 



5-2p 



dp 



+ Kb 



-6 + 2p 



d_ 
dp. 



(6.12) 
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which is the same result we obtained in eq. (|6.4|). Eqs. (|6.10| )- (|6.12| ) agree with eqs. (|5T^)- 
( p.7| ) on setting Gb — > 0. One can also understand the origin of the 5-function terms in 
sec. 5. Since the decay distribution does not vanish at the end point, the derivatives in 
eq. ( |6.11| ) produce 5-functions and derivatives of 5-functions. 

One can also apply the averaging method to obtain the spin-dependent form-factors. 
One considers the b-quark in the hadron to have a distribution of spin and momentum, 
with 

{Sn = il + eb + Kb)s^, 

{S^k") = TUbKbS^v" + IrubKbS^v^', (6.13) 
{S^k'^k^) = -'^mlKbig"'^ - v"vf^)s^' 
where S'^ is the quark spin, and is the hadron spin. 



7. Decay Distributions for Hadrons Containing a c-Quark 

The decay distributions for semileptonic c-quark decay can be readily obtained from 
the calculations in the previous section for b-quark decay. The charged lepton distribution 
in c-decay is equal to the neutrino distribution in b-decay, and vice-versa. This is equivalent 
to changing the signs of Ws, 6*3, Gs and Gg. The lepton energy spectra are less singular 
than for b-decay, because the free quark decay rate vanishes at the endpoint. The double 
differential c-decay distribution for an unpolarized hadron He containing a c-quark is 

e{z) \ 12y{l -p-y) + l2yE,{l + p) - 2AyK, + ^yG,{p - 2] 



Vc dy dq^ 



y I (7.1) 

+ 4K,(3^4 - 6fy - AqS + ^fy^ + y^) 

+ 8g2Ge(l - y){e + ?/ - + 5'{z)Kj^{y^ - ffil - y), 

J y 

where 

r.Hv-,/^. (7.2, 
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Integrating this with respect to and then y gives 



1 dV 



12(l-y)y2_2Vp + 



2^2 



I2y'p 



Tc dy \ (1 - y) 

+ EAl2{2-y)y'+^'^" 



2 , 12y2(2/-2)p2 



(1-2/)^ 
8y^(3-2y)p^ 



3 , 82/3p2 



and 



1 - 80 + 



.3 A 



p*-12p^logp 



+ £;J5-24p + 2V 



+ 3p^-12p2 logp 



+ + 32p - 24p^ - 2p* + 24p^ log p 



+ GJ-2 + 16p - 16p3 + 2p4 + 24p2 bgp I 



The p ^ hmits of eqs. — ( |7^ are 



1 rfr 



^(^) 12(2; - 2/2) + 122/S, - 24yK, - 16yG, 



1 



+ 5(z)-<^ 12£;er(2/ - l)(r - 22/ + 

y 



- AK,{-?,ct + Wy + Wy - QqY - yn 

fG^iy -l){f + y- /)) + 5'{z)K^%{y' - ffil 



j_dr 

Tc dy 



12{l-y)yH+EJl2{2-y)y 



+ KJ -82/2(3 + y)+ 45(1 -y)\+Gj -8y 



T = To{l + 5E,-6K,-2Gc} 
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The cos 9 terms in the decay distributions for a polarized Ag are 



- 6y(l - p)Ke} + 5{z)^2K,{-t - Sfy + 3qY + V^) 



1 dr 

Fc dy dcos 9 



y 



+ 



cos 9, 



(l + 6e)^6(l-y)y2-12yV+^^ 



+ + (1-^)3 



cos 0, 



1 

Ff, cZcos ^ 



= . . . + (1 + ee - Kc) 



The p — > hmits of these distributions are 



1 dT 



. + 



9{z) (1 + ee) 62/(1-2/) 



Fc dy dq^ dcos 9 

- 6yK,] + d{z)l2Ka{-q'' - Sfy + 3qV + 2/') 

cos 9, 



+ 5'{z)K,%{y'-ef{l-y) 



y 



1 dT 
Fc dy dcos 9 



. + 



(l + ee){6(l-2/)2/'} 



+ Ke{-102/' + 25(l-2/)} 



cos 6*, 



1 dV 
Fc dcos ^ 



...+ -(l + ec-i^c)cos^. 



(7.8) 



(7.9) 



cos 6*. (7.10) 



(7.11) 



(7.12) 



(7.13) 



8. Inclusive H^j Xy^ev^ near the boundary of phase space 

For an exclusive decay Hf, — > XeUf, the kinematically allowed region of phase space is 



0<q^ < 2E^Mh, + 



{2Ee-MH,)' 



(8.1) 
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The maximum value of the electron energy E^^^^^ occurs when the right hand side of 
eq. is zero, 

, , Ml -Ml 

The results of sec. 5 show that in QCD, the inclusive decay kinematics are governed by 
the quark mass mf,, rather than the hadron mass Mh^,- The kinematically allowed region 
in the Dalitz plot is 

Q<q^ < 2E,m^ + -— ^ 8.3 

{2E^ - mi,) 

where rrij is the charm quark mass for b —>■ c transitions and zero for h ^ u transitions 
(neglecting light quark masses). The difference between these two kinematic regions is 
shown in fig. 4 for h ^ u decay. 

For h ^ u transitions the region (8^) of the q'^^Ee plane becomes the interior of 
a right triangle with sides = 0,q^ = 2Ef,mh and E^ = rnb/2. Since is less than 
Mh^ a comparison of eqs. ( |8.3| ) and ( |8.1| ) reveals that a part of the kinematically allowed 
phase space is not populated. This part corresponds to the production of states with mass 
squared 

ML ^ (Mh, - m,)MH, (^1 - ^ j . (8.4) 

The physical origin of this discrepancy has been discussed by Isgur, Scora, Grinstein and 
Wise and by Isgur [|1^. For simplicity, consider QCD in the large Nc limit where the 
final state X must be a qq bound state (nonresonant final states are produced with an 
amplitude suppressed by a factor of l/v^iVc). Provided /Mfj^ is not very close to unity, 
the production of states with masses much less than in eq. 0.4| ) is strongly suppressed 
by hadronic form factors. For very large rrih this form factor suppression arises from the 
transfer of a large momentum to the spectator antiquark by a single hard gluon. Clearly 
such effects are subdominant to those given by eq. ( |5.8| ). However, when the mass of 
the final qq resonance becomes of order a/^^t^Aqcd (see eq. ( |8.4D ), the antiquark in the 
X meson has such a broad distribution of momentum that no form factor suppression is 
required. Thus the differential decay rate in eq. (|5.8D corresponds to a sum over exclusive 
final states with masses greater than ~ lym^jAqcD- Note that for large mi, these hadronic 
masses are much greater than the QCD scale. 

There are theoretical limitations on the extent to which our prediction for the differ- 
ential Hi, — * XueVe decay rate dV / dq^dEf, can be compared with experiment. In the region 
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of phase space very near = Mfj^, which corresponds to low mass final hadronic states 
recoiling at low momentum, our expression for the differential decay rate is not valid (the 
operator product expansion cannot be justified for low mass states). Also, the appearance 
in eq. ( |5.8|) of a delta function and its first derivative indicates that along the boundaries 
= 2mi,Ee and E^, = m{,/2 the differential decay rate must be smeared over a region 
of electron energies. The amount of smearing necessary is determined by demanding that 
corrections proportional to Kh,Ei, and Gi, give a contribution to the smeared differential 
rate that is small compared with the leading "free quark decay" contribution. 

For definiteness consider the boundary = 2mi,Ef,. In the free b — ueVe decay it 
corresponds to a configuration (in the b rest frame) where the electron and u quark are 
moving in the same direction and the antineutrino goes in the opposite direction. Since 
the weak current is left handed the free quark decay amplitude vanishes here by angular 
momentum conservation. Define (recall y = 2Ef,/mi,, = q^ /ml) 

S{q') = J dy^^^Wiy), (8.5) 

where 

1 e-(^-('^'+^))'/^'. (8.6) 

S{q'^) corresponds to smearing the differential cross section at y = q'^ + e over a region 
of y of order e. (The averaging is aX y = q^ + e instead of y = q^ so that the term 
containing a derivative of a delta function is not made artificially small.) Demanding 
that the leading free quark decay contribution to S{q'^) be large compared to that of the 
corrections proportional to Kb, Et and Gb gives the condition 

> (AQCD/m,)2. (8.7) 

This corresponds to smearing over a range of electron energies AEe > Aqcd- In deriving 
eq. ( |8.7] ) the matrix elements Kb, Eb and Gb were estimated to be of order (AgcD/^b)^- 
A similar result holds for the region near the boundary of phase space at E^ = mb/2. The 
differential cross section in eq. (|5.8|) must be smeared over a region of electron energies 
AEe > Aqcd to be physically meaningful. Away from the boundaries Ee = mb/2 and 
q^ = 2mbEe of phase space some smearing may also be required. As one varies E^ (at fixed 
q'^) new thresholds are encountered corresponding to different Mx in eq. (|8.1| ). However, 
here the situation is likely to be similar to R{s) = a{e~^e~ hadrons)/cr(e"'"e~ fJ.~^fJ.~) 
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where at large s the density of hadronic states is so large that the size of the region s that 
must be averaged over becomes negligibly small. 

In semileptonic B decay the endpoint region of electron energies, where B XcCVf, 
decays are forbidden, is of great interest because of its potential use to determine \Vub\- 
However, in this endpoint region 2.2 GeV < E^, < 2.6 GeV our predictions are not very 
useful. Very near = our results don't apply and away from this value of the 
averaging of electron energies over a region Ai?e > Aqcd makes it impossible to isolate 
the endpoint part of the electron spectrum. This is discussed more quantitatively in the 
next section. 

9. Numerical Results 

The corrections to the decay distribution for semileptonic Hi, decay have been 

computed in terms of the parameters Gh and Ki, . These parameters depend on the hadron 
ifb, and will be denoted by Gi,{Hh) and Ki,{Hi,). The parameter Gb is the leading operator 
that breaks the heavy quark spin symmetry, so its matrix element can be determined 
in terms of the hyperfine spin-splittings within heavy quark multiplets. For example, 
GbiAh) = and 

rribGbiB) = [M{B*) - M{B)] , 771,^,(5*) = ^ [M{B*) - M{B)] . (9.1) 

The experimental value of 46 MeV for the B* — B mass difference determines Gb{B) 
to be —0.0065 (where nib can be equated with the hadron mass to this order). The 
parameter Kh{B) cannot be simply determined. Quark model estimates suggest that 
Kb{B) is approximately 0.01. The size of the l/mf correction to the lepton spectrum 
dV/dy for B X^eV^ decay is plotted in fig. 5 for Gb{B) = -0.0065 and KbiB) = 0.01. 
The plot shows the ratio of the distribution eq. ( ^.6|) to the free quark decay spectrum 
without any smearing. Fig. 6 shows the same result as a percentage correction, so that 
one can see how the corrections become large only near the endpoint. The results must 
be averaged in y over a region sufficiently large that the correction terms are small. The 
1/m^ corrections to charm decay are a factor of (m^/mc)^ larger than the corrections for 
b-decay. 

To better understand the size of the 1 / corrections, it is interesting to plot the B —>■ 
XcCVf, lepton spectrum including 1 / corrections without smearing. This is illustrated in 
fig. 7, where the free quark b ^ c decay and h u decay electron spectra are compared to 
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the B — s> XcCVe decay spectrum including 1/m^ corrections. The peak near the endpoint 
for B — > XcCVe decay gives an indication of the minimum size of the smearing region that 
must be used before the QCD calculation is valid. The peak must be smeared over a large 
enough region that it produces a small correction to the decay spectrum. This indicates 
that smearing region should be at least Ay = 0.2, which corresponds to a lepton energy of 
around 500 MeV. 

The dominant uncertainty in the extraction of Vub is the shape of the endpoint of 
the B —>■ XuCVe lepton spectrum. It is difficult to estimate the size of the smearing 
necessary in QCD from the unsmeared spectrum eq. ( |5.10| ) since it contains 5-function 
singularities at the endpoint. Figure 8 shows the smeared lepton decay distribution using 
Gaussian smearing in y with different widths. The smearing extends the spectrum beyond 
the parton model endpoint y = 1, but the curves are only plotted for < y < 1. For 
simplicity, only the 5-functions in eq. ( |5.10| ) have been smeared, and then added to the 
remaining terms. Clearly, the dip in the spectrum for a smearing width of Ay = 0.1 is 
unphysical. The curves in fig. 8 indicate that the minimum smearing width in Ay is around 
0.2, which corresponds to a lepton energy width of 500 MeV. 

Even though the parameters Kb[Hb) are not known except by model calculations, one 
can determine the difference of between two hadrons. For example, the masses of the 
Kb baryon and the B and B* mesons are 

M{Kb) =mb + A(A6) + mbKb{Kb). 

M{B) = m6 + A(S) + mbKb{B) + mbGb{B), (9.2) 
M{B*) =mb + A{B) + rubKbiB) - ^mbGb{B), 

including all corrections to order l/rub in the mass. We have used the heavy quark spin 
symmetry relations Gb{B) = —3Gb{B*) and Kb{B) = Kb{B*). The difference between 
the Ab mass and the average meson mass M{B)avg = {M{B) + 3M{B*))/4 is 

M{Ab) - M(S)avg = A(Ab) - A(S) + {Kb{Ab) - Kb{B)) . (9.3) 

A similar expression holds for the D mesons and the A^ baryon. Heavy quark symmetry 
implies that A{Hc) = A{Hb), and mlKb{Hb) = mlKc{Hc), so that 

[M(A,) - M(I^)avg] - [M{Ab) - M{B),,^] =mb(—- l) {Kb{Ab) - Kb{B)) . (9.4) 

33 



Equating the quark masses with the meson masses to this order gives Ki,[A.i,) — Ki,{B) = 
—0.002 ± 0.006, using the present value of 5641 ± 50 MeV for the A;, mass. 

The differences between Gb and Kb for different hadrons gives predictions for the 
differences in the semileptonic decay distributions. The order as corrections due to gluon 
radiation that have been computed [000 are the same for aU hadrons, since they correct 
the free quark decay formula, and cancel in the difference. Thus the difference is known to 
corrections of order as(m&)/m& and For example, one can compute the differences 

in the total semileptonic decay widths for Hb XueV^ decay, 

riff!) + r(g^) = I ^^'^"'^ " ^'^"'^^ " I ^^'^^'^ " ^'^^'^^ ■ ^^-^^ 

This gives 

T(Ab)-T(B) 9 (M(B*) - M(B)) , , , , , 

^ ' = 4 - (^^(^'') - ^''(^)) = 0-018 ± 0.006. (9.6) 

Similarly, the decay width differences for He —>■ XdC^Ue decay are 
r(Ac)-r(D) 9 {M{D*) - M{D)) 



r{D) 4 M{D) 



{K,{A,)-K,{D))=0.16±0M. (9.7) 



The uncertainties in eqs. ( |9.6| ) and ( |9.7| ) are due to the ±50 MeV uncertainty in the Af, mass. 
The differences of the decay widths for the decay modes Hb — > X^eVe and H^ Xge'^Uf. 
can be computed similarly using eqs. ( |5.11| ) and ([7.4|). 



10. Conclusions 

One of the important points of the analysis in this paper is that the decay distribution 
is determined by the quark mass m^, rather than the hadron mass Mhi,- Thus the free 
quark decay rate depends on m^, not Mfj^. The difference between the two masses 
is the A parameter of the heavy quark theory. There should have been l/rub corrections 
proportional to A if a free quark decay model with the decay rate given by Mh^ was 
appropriate. Corrections of this form are absent in QCD. 

The endpoint region of the inclusive lepton spectrum in semileptonic b-decay is impor- 
tant for the extraction of V^b- The QCD calculation near the endpoint must be smeared 
over a large region (around 500 MeV) before it can be compared to experiment. This 
region is larger than the difference in the endpoints of the b u and b c decays. This 
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means that the extraction of V^b stiU requires modehng the endpoint. The QCD computa- 
tion does provide some constraints on the model. Any model which has l/m^ corrections, 
or has a decay distribution given by hadron kinematics instead of quark kinematics, is in 
contradiction with QCD. 

We have presented results for decays of hadrons containing a b or c quark. Formally, 
perturbative corrections are of the form as{rnb) or asirric)- However, since the final state 
in free quark decay is three-body, perturbative corrections to the total semileptonic decay 
rate may be better represented by as(mb/3) or as(mc/3). Only a higher order perturbative 
calculation can resolve this issue. Because of this, we don't have confidence in applying 
the results of sec. 7 to the decays of hadrons containing a c-quark. 
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